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Abstract. Let U he a matrix chosen randomly, with respect to 
Haar measure, from the unitary group U{d). We express the mo- 
ments of the trace of any submatrix of {/ as a sum over parti- 
tions whose terms count certain standard and semistandard Young 
tableaux. Using this combinatorial interpretation, we obtain a 
simple closed form for the moments of an individual entry of a 
random unitary matrix and use this to deduce that the rescaled 
entries converge in moments to standard complex Gaussian ran- 
dom variables. In addition, we recover a well-known theorem of 
E. Rains which shows that the moments of the trace of a random 
unitary matrix enumerate permutations with restricted increasing 
subsequence length. 



1. Introduction 

Consider the unitary group U{d) as a probability space under nor- 
malized Haar measure. Given a random variable X : U{d) C, its 
expected value is defined to be 

Ec/(,)(X) = / XdU. 
Ju{d) 

When studying a random variable X, one often wishes to know its 
moments 

since in many situations the moments of X uniquely determine its 
distribution. X is called a polynomial random variable if there is a 
polynomial / G C[a;ii, . . . , Xdd\ such that 

X{U) = /(nil, • • ■,Udd) 

for all U G U{d). If X is a polynomial random variable, then 

Et;(rf)(X'"X") =0 if m^^n, 
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(see [Ij), so knowledge of the moments of X reduces to knowledge of 
the quantities 

Recently, certain polynomial random variables on the unitary group 
have been shown to possess interesting combinatorial properties. For 
a random matrix U G U{d), write its characteristic polynomial as 

d 

det{U-zI) = {-irY.{-iySc,{U)z'-^. 

j=0 

Scj{U) is called the j-th secular coefficient oi U. In particular, 

Sci(t/) = Tr([/) 

Scd([/) = det([/). 
It was shown by Rains in [5j that 

E^w(|Tr(f/)P")= 

AI-n,£(A)<d 

the number of pairs of semistandard Young tableaux on the same Young 
diagram, over all possible Young diagrams with n boxes and at most 
d rows. By the Schensted correspondence, this is equal to the number 
of permutations in the symmetric group Sn with no increasing subse- 
quence of length greater than d. 

The other secular coefficients of a random unitary matrix also en- 
code interesting combinatorial information, as shown by Diaconis and 
Gamburd in ^2]. For jn < d, [2j proves that 

Euid){\Sc,{U)n=H4j), 

where Hn{j) is the number of n x n matrices whose entries are non- 
negative integers and whose rows and columns all sum to j ("magic 
squares"). The method of proof given in |JJ is simple and elegant, and 
consists of two main ingredients: 

• The secular coefficients of U are the elementary symmetric func- 
tions apphed to the eigenvalues of U. 

• The elementary symmetric functions can be written as linear 
combinations of Schur functions, which are the irreducible char- 
acters of the unitary group and thus satisfy orthogonality rela- 
tions. 

Now suppose that we are given a random matrix U & U{d), and we 
want to calculate the moments of a single entry of U 
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(note that since the permutation matrices are in U{d) and Haar mea- 
sure is translation invariant, all entries of U are equidistributed). More 
generally, for a positive integer k with 1 < k < d, let Uk denoted the 
k X k upper left corner of U. We could ask about the moments of the 
matrix Uk 

E^w(|Tr(f/,)r") 

(again, since the permutation matrices are in U{d), the traces of any 
two k X k submatrices of U are equidistributed). The methods used in 
jS], PI will not work in this situation, since we are no longer dealing 
with the eigenvalues of a unitary matrix, but rather the eigenvalues of 
submatrices of a random unitary matrix. However, quite surprisingly, 
there is a simple combinatorial formula for these moments. 

Theorem 1. Let U be a matrix chosen randomly with respect to Haar 
measure from the unitary group U{d), and let Uk be its k x k upper left 
corner. We have 

Eu^,)imuk)n= y. if'^^''^'^^^ 



SXdil) 

Ahn, e{X)<k ■ ^ ' 

Here, for a positive integer r, SA,r(l) denotes the number of semistan- 
dard Young tableaux on the shape A with entries from the set [r] = 
{l,...,r}. 

The notation SA,r(l) is shorthand for the value of the Schur function 
s\ obtained by setting the first r variables equal to 1 and making the 
remaining variables 

sa(1,1,.. .,1,0,0,...). 

The matrix f/^ is called a truncation of U. Truncations of random 
unitary matrices where first studied in [7j from an analytical point 
of view, where it was shown that Uk is a contraction (i.e. all of its 
eigenvalues lie in the closed unit disc in C), and the joint probability 
density function of the eigenvalue was found to be 

i<j j=l 

Cd,k is a normalization constant that was found in T to be 

j=o \ J y 

Note that Rains's result is an immediate corollary of Theorem 1; it 
is simply the case k = d. 
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Corollary 1.1. In the special case k = d, we have 

Ahn,£(A)<d 



Proof. When k = d, 



SA,d(l) SA,d(l) 



□ 



Rains's theorem is the extreme case k = d oi Theorem 1. The 
other extreme = 1 is also of interest, since this corresponds to the 
computation of the moments of a single entry of a random unitary 
matrix. 

Corollary 1.2. In the special case k = 1, we have 

^ /, ,9r,N fd + n — l\'^ 

^' d{d+l) . . .{d + n-l) \ n J 

Proof. For k = 1 the only contribution to the sum is made by the single 
partition whose diagram is a row of n boxes. Thus we have 



lu{d) 



The generalized hook length formula asserts that 

□6A ^ ^ 

where c(n) is the content of the box, and h{\Ji) is its hook length (see 
jH]). For the single row partition of n, this gives 

d{d+l)...{d + n-l) 
= , 

and the result follows. □ 



Explicitly knowing the moments of Uij makes it easy to determine its 
limiting distribution. Recall that if x, y are Gaussian random variables 
with mean and variance 1/2, then the random variable z = x + iy is 
called a standard complex Gaussian. 

Corollary 1.3. As d ^ oo, the random variable Vduij converges in 
distribution to a standard complex Gaussian random variable. 



TRUNCATIONS OF A RANDOM UNITARY MATRIX AND YOUNG TABLEAU}© 



Proof. It is well-known (see for instance jlj) that the moments of a 
standard complex Gaussian z are given by 

Corollary 1.2 shows that 



= 6.. 



d{d+ 1) ...{d + n-1) 



□ 



In order to prove Theorem 1, one needs to connect unitary expec- 
tations to symmetric function theory by some method other than ap- 
plying symmetric functions to eigenvalues. This can be done using 
the Weingarten function introduced in jl], which is a powerful tool for 
computing the moments of polynomial random variables on the unitary 
group. The Weingarten function has already been used in free proba- 
bility theory to prove asymptotic freeness results for random unitary 
matrices (see the recent book jS] for a clear account of this). 



2. The Weingarten Fungtion 

For any positive integers d and n, define a function Wg{d, n, ■) : Sn 
Qby 

Ahn, e{\)<d ^'"^^ ' 

where is the irreducible character of Sn labeled by A. 
The following integration formula was proved in P . 

Theorem 2. Let i, , j' : [n] — ^ [d] be any functions. Then 

■ ■ ■ Ui(^n)j{n)Ui(l)j(l) ■ ■ .Ui(n)j(n)dU 

U{d) 

= E • • • ^iin)i'(ain))Sj(l)j'(r(l)) ■ ■ ■ Sj(n)j'(T(n)) Wg((i, U, Ta~^), 

where 6 is the Kronecker delta. 

We can succinctly express the moments of Tr(f/fc) in terms of the 
Weingarten function as follows: 
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Lemma 1. For any positive integers n and k, where 1 < k < d, we 
have 

where the outer sum runs over all weak k-part compositions a ofn, and 
the inner sum runs over all permutations in the Young subgroup Sa of 

Sn- 

Proof. This is really just a calculation. We expand 



Tr{Uk)\^^=\un + --- + Ukk\^^ 

= {Uu H h H h UkkY 



a /3 

where we are summing over all pairs of weak fc-part compositions 
of n, a = {ai, . . . , a^) and (3 — {bi, . . . ,bk). We are using multi- index 
notation, 



n\ / n 



aj \ai,...,ak 
u^^ut\...ull 

n\ f n 
(5) ^ W...,bk 



Hence 



We will use the Weingarten integration formula to evaluate the ex- 
pectation Ku(^d){U'°'u^) for a fixed pair of compositions a, (3. Implicitly 
define coordinate functions ia, ja,i/3, jp '■ [n] — > [d] by setting 

/ Ui^{l)ja{l) ■ ■ ■ Uic(n)jc{n)Uig(l)j0(l) . . . Ui^{n)ji3{n)dU 
JU{d) 

Ju{d) 
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= ^ Si^{l)i0{a{l)) ■ ■ ■ Si^{n)if3{ain))Sj„{l)j0{T{l)) ■ ■ ■ (^J„(n)j^(r(n)) Wg(d, n, TO"-^). 

Since we are only taking entries from the diagonal, we have ia = ja 
and ifs = ji^. Moreover, the level sets of these functions are easy to read 
off: 

= [ai + 1,01 + 02] 
C^(^) = [n-ak + l,n] 

= [n-hk + l,n] 



and 



Hence 



ai ai 

r=l r=l 

ai+a2 ai+a2 

r=ai+l r=ai+l 



r=n— a^+l r=n— a^+l 

Thus in order for the product 

^ic{l)i0{a{l)) ■ ■ ■ (^ia(n)i^((T(n)) 

to be nonzero, we see that a must bijectively map the interval [1, oi] 
onto the interval and a must also bijectively map the interval 

[oi + 1, oi + 02] onto the interval [61 + 1, 61 + 62], etc. Similarly, in order 
for the product 
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to be nonzero, wc sec that r must bijectively map the interval [1, ai] 
onto the interval [l,&i], and a must also bijectively map the interval 
[oi + 1, ai + 02] onto the interval [61 + 1, foi + 62], etc. Thus we see that 
the expectation Ku(^d){u'^u^) is zero unless: 

• a — P, i.e. these two are the same weak k-part composition of 
n, 

• a,T are both in the Young subgroup Sa, i.e. the subgroup of 
permutations in Sn that permute the first Oi symbols amongst 
themselves, the next 02 symbols amongst themselves, etc. 



Thus 



a,T€Sa 

2 



□ 



3. Proof of the Main Theorem 
We are now in a position to prove Theorem 1. 
Proof. Having proved that 

Eu^a){\Tr(U,)n ^ n\J2 E Wg(d,n,a), 



(JGSc, 



we will work with the sum on the right. Plugging in the definition of 
the Weingarten function, we have 

"'E(:)Ewgw»..)^»!E(:)E;^ E St)^V). 

Changing order of summation, this becomes 

E (/-v^j^E^Ex^w- E (/^)^j^E(i.xV. 
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where the inner product {■,-)sa is the averaged dot product on the 
space CF{Sa) of complex-valued class functions on the group Sa- Note 
that this sum may be written as 



Ahn, ii\)<d ^^•'^'^^^ . 

where x'^ isZ restriction of the irreducible character x'^ of Sn to 

the subgroup Sa- Now, the function which is identically 1 is the char- 
acter of the trivial representation of Sa- Thus we may apply Frobenius 
reciprocity (see for instance P); 

where 1 is the induction of the trivial character of 5*^ to Sn- 

The final step in the proof relies on the characteristic map ch" : 
CF{Sn) —>■ A", where A" is the inner product space of degree n sym- 
metric functions equipped with the Hall inner product (■, ■)a" (see for 
jH] or jH]). The class function is an isometry, and has the following 
important properties: 

Ch„(l T£) = ha 

ch„(x^) = Sa, 

where ha is the complete homogeneous symmetric function indexed by 
a, and s\ is the Schur function indexed by A. Thus we have 

(1 Tf:,X^)s„ = (ch„(l ll:),cK{x^))s„ = {ha,Sx)Ar^- 

It is well-known that the Schur functions constitute an orthonormal 
basis for A", and that the coordinates of the complete homogeneous 
symmetric functions with respect to the basis of Schur functions are 
the Kostka numbers (see 0). That is, 

ha ^ ^ K ^-1 

fi\—n 

where the Kostka number K^j^a is by definition the number of semis- 
tandard Young tableaux on the diagram of /i with content vector a. 
Thus, 

{ha,Sx)A^^ = y^^K^a{Sf_i,Sx)A^ = Kx,a- 
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Thus we have 



Xhn, i{\)<d I \ / a Ahn, E{X)<d > \ / q. 



gA,fc(l) 

Ahn, £(A)<d ^^■"'•^^^ 



where the last equahty follows from the fact that, by definition, 

"^Kxa = Sa,/c(1), 

a 

since the sum runs over all weak fc-part compositions of n. 

Finally, we remark that if A is a partition of n with i[X) > k, then 
SA,fc(l) = 0. Hence, 



>S\,k{l) _ ^^AN2'SA,fc(l 

which proves our theorem. 



Ahn, i{X)<d ^^''^^ ' Ahn, £(A)<fc ^^''^^^ 



□ 



4. Conclusion 

In this paper, we have only investigated the moments of the trace of 
a truncation of a random unitary matrix. It seems possible to analyze 
the moments of the other secular coefficients of a truncation by the 
same method, and it would be interesting to see what combinatorial 
interpretations can be given for the moments of these coefficients. 

In |2] , the moments of secular coefficients of random orthogonal and 
symplectic matrices are investigated, and results analogous to the uni- 
tary case are proved. In [J, a notion of Weingarten function is defined 
for the orthogonal and symplectic groups. It is likely possible to analyze 
the moments of secular coefficients of truncations of random orthog- 
onal and symplectic matrices using the Weingarten function for these 
groups. 
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